Kondo Shuttling in Nanoelectromechanical Single-Electron Transistor 
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We investigate theoretically a mechanically assisted Kondo effect and electric charge shuttling in 
nanoelectromechanical single-electron transistor (NEM-SET). It is shown that the mechanical mo- 
tion of the central island (a small metallic particle) with the spin results in the time dependent tun- 
neling width r(i) which leads to effective increase of the Kondo temperature. The time-dependent 
oscillating Kondo temperature 2V(i) changes the scaling behavior of the differential conductance re- 
sulting in the suppression of transport in a strong coupling- and its enhancement in a weak coupling 
regimes. The conditions for fine-tuning of the Abrikosov-Suhl resonance and possible experimental 
realization of the Kondo shuttling are discussed. 

PACS numbers: 73.23.Hk, 72.15.Qm, 73.63.-b, 63.22.-|-m 



Kondo resonance tunneling predicted in ll| and ob- 
served experimentally [3| attracts a great deal of cur- 
rent attention as a possible base for manipulating spin 
transport. Nanomechanical shuttling (NMS) [2| offers a 
unique platform for design of a single electron transistor 
(SET) in which spin switch/transfer can be controlled 
electromechanically; the first successful experimental im- 
plementation of the NMS were reported in |j, |5j . In this 
Letter we develop a theory of a nanomechanical shuttling 
device that utilizes Kondo resonance effect (KR) and thus 
breaks the ground for a new class of effects integrating 
both phenomena. 

The Kondo effect in electron tunneling results from the 
spin exchange between electrons in the leads and the is- 
land (quantum dot) that couples the leads and manifests 
itself as a sharp zero bias anomaly in low temperatures 
tunneling conductance. Many-particle interactions and 
the tunneling renormalize the electron spectrum enabling 
KR both for odd [2| and even laQ electron occupations. 
In the latter case the KR is caused by the singlet-triplet 
crossover in the ground state (see |8J for a review). 

A general shuttle mechanism for a charge transfer de- 
scribed in [3J implies periodic charging and de-charging 
of the oscillating nano-particle. As the bias exceeds some 
threshold value, the shuttling particle oscillates with the 
constant amplitude along a classical trajectory. A model 
for a nano-electromechanical single electron transistor 
(NEM-SET), where either a small (nano-scale) metallic 
particle or a single molecule oscillate under the external 
time-dependent electric field was studied theoretically in 
|9j-[l^. In experimental realizations, the moving dot was 
mounted as a electromechanical pendulum formed by a 
gold clapper p, or a silicon nanopillar |j|. Single elec- 
tron tunneling in a molecular conductor with the center 
of mass motion was realized in |l5l Il6| . The vibration 
induced Kondo effect in metal-organic complexes was fur- 
ther explored in [Tg ■ Experimental realization of Kondo 



effect in electron shuttling is being one of the most chal- 
lenging tasks of current nano- physics. 

In this Letter we investigate the effect of the spin de- 
grees of freedom on the single electron transport through 
the NEM-SET. We show that mechanical shuttling of 
a nano-particle between the leads causes the sequential 
in time reconstruction of its tunnel electronic states and 
gives rise to Kondo effect. By the analogy with the effect 
of the sequential in time recharging of the particle in or- 
dinary charge shuttling effect [3| we call this phenomenon 
Kondo shuttling. 

Building on the analogy with shuttling experiments of 
1^ |5J J we consider the device where an isolated nanoma- 
chined island oscillates between two electrodes. We, how- 
ever, are interested in a regime where the applied voltage 
is low enough so that the field emission of many electrons, 
which was the main mechanism of tunneling in those ex- 
periments, should be neglected. Note further that the 
characteristic de Broglie wave length associated with the 
dot should be much shorter than typical displacements 
allowing thus for a classical treatment of the mechani- 
cal motion of the nano-particle. The condition MI^Tk, 
necessary to eliminate decoherence effects, requires for 
e.g. planar quantum dots with the Kondo temperature 
Tfj-^lOOmii', the condition Q,<\GHz for oscillation fre- 
quencies to hold; this frequency range is experimentally 
feasible [j, [Sj- The shuttling island then is to be con- 
sidered as a " mobile quantum impurity" , and transport 
experiments will detect the influence of mechanical mo- 
tion on a differential conductance. If the dot is small 
enough, then the Coulomb blockade guarantees the single 
electron tunneling or cotunncling regime, which is neces- 
sary for realization of Kondo effect [ll Q • Cotunneling 
process is accompanied by the change of spin projection 
in the process of charging/discharging of the shuttle and 
therefore is closely related to the spin/charge pumping 
problem |l8l |. 
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FIG. 1: a) Nanomechanical resonator with the odd number 
of electrons as a "mobile quantum impurity", b) Time evolu- 
tion of the level spacing for the shuttle with the even number 
of electrons measured with respect to evolving Kondo tem- 
perature, c) Pulse ZBA in tunneling conduction for Kondo 
shuttling with the singlet/triplet transition. 



We apply our theory to planar quantum dots in semi- 
conductor heterostructures [l^. In these systems the 
KR tunneling may be realized both for odd and even 
electron occupation A^ |3, la| ■ If iV is odd, the last oc- 
cupied level in the island is occupied by single electron 
(Fig.l.a). In this case the Kondo effect occurs and the 
corresponding Kondo temperature, T]^, is that of a static 
regime. The question we wish to address is how Tr- is 
influenced by the adiabatic mechanical motion of the is- 
land. If, on the contrary, TV is even, the ground state 
of the nano-particle is singlet S'=0 whereas the excited 
state is triplet S—1 (Fig.l.b). Under the assumption 
that the energy difference between the singlet and triplet 
states AsT=S—Jex exceeds the actual Kondo tempera- 
ture Ast>Tk, the Kondo effect is absent in the static 
limit (here (5=62— ei and Jex is ferromagnetic exchange, 
ei=i,2 are the energies of single-electron states for the 
last occupied and first empty levels). However the me- 
chanically induced Kondo effect may arise provided the 
inequality AsT<TKit) is satisfied due to mechanical mo- 
tion of the dot. 

Now we turn to a quantitative description of Kondo 
shuttling. The Hamiltonian H=HQ+Htun is given by 

Hq = '^Ska^acl^^a^ka.a +^['^i ~ e£x]dl^dicr + Un^ 

k.a ia 

Htun = Y. T^''^^M.,c.d^a + H.C], (1) 

ika.a 

where c\^^ d\^ create an electron in the lead a~L,R, or 

the dot level £4=1,2, respectively, n—^^^dj^dia, £ is the 
electric field between the leads. The tunnelling matrix el- 



ement Tl^'^{x)= T£ '^^exp[=Fa;(i)/Ao], depends exponen- 
tially on the ratio of the time-dependent displacement 
x{t) (which is considered to be a given harmonic func- 
tion of the time) and the electronic tunnelling length Aq. 
We begin with the discussion of an odd N, S = 1/2, 
case (Fig.l.a). Then only the state with i=l retains in 
(^, and hereafter we omit this index. In order to find 
an analytic solution, we assume that if x{t) varies adi- 
abatically slow (on the scale of the tunneling recharg- 
ing time), there is no charge shuttling due to multiple 
recharging processes Q, but the KR cotunneling oc- 
curs. The time -dep endent tunneling width is Ta(t) = 
2?: pQ\Ta{x{t))\'^ [20|, where po is the density of states at 
the leads Fermi level. The adiabaticity condition reads: 
=min[Vri(ir 



UVKTr <r, with r=min[^r2(i)+r|j(i)]. We apply 
the time-dependent Schrieffer- Wolff transformation and 
obtain the time-dependent Kondo Hamiltonian [20j as 

H = Hq+ ^ Jaa'{t)[ffaa'S + -S^„>]cl^^^Ck'a',a' 
kaa.k' a' a' 

, (2) 

where Ja,a'(i)=yr40lVW/(7rpo^d(t)) and 

S=\d\,Saa'dcr' ■ Without the loss of the generality, 
we can restrict ourselves to symmetric Anderson model 
{eo~—U/2 in the static limit) and neglect the renormal- 
ization of a single electron level position by tunneling 
(r/[/^l) and also by its shift from the equilibrium 
position [(ef ^)/C/<l] such as Ed{t)fvE^=U/A. 

As long as the nano-particle does not subject to the ex- 
ternal time-dependent electric field, the Kondo tempera- 
ture is given by T]^=£)oexp[— (7rC/)/(8ro)] (for simplicity 
we assumed that rL(0)=rfl.(0)=ro; i:'o=\/2roC//7r plays 
the role of effective bandwidth). As the nano-particle 
moves adiabatically, ftri<cr, the decoherence effects are 
small provided ftfi^T^ (see the discussion below). In 
this case the time can be treated as an external param- 
eter, and the renormalization group equations for the 
Hamiltonian can be solved the same manner as those 
for the equilibrium |20j. As a result, the Kondo temper- 
ature becomes time oscillating: 



TK{t) = D{t)c^p 



■kU 



8roCOsh(2a;(t)/Ao) 



(3) 



Neglecting the weak time-dependence of the effective 
bandwidth D{t)^DQ, we arrive at the following expres- 
sion for the time-averaged Kondo temperature: 



(Tk)=TO( exp 



ttU sinh^(2;(t)/Ao) 



4rol-h2sinh^(a;(i)/Ao) 



(4) 



Here (...) denotes averaging over the period of the me- 
chanical oscillation. The expression Q acquires espe- 
cially transparent form when the amplitude of mechanical 
vibrations A is small: A<Xo. In this case the Kondo tem- 
perature can be written as (T/^ )=T]^exp(— 21^), with the 
Debye-Waller-like exponent VF=-7rC/(a;2(t)))/(8roAg), 



giving rise to the enhancement of the static Kondo 
temperature. This counter-intuitive exponentially large 
Debye- Waller factor results from the strong asymmetry 
in the tunneling rate at the turning points of the nano- 
particle trajectory. Note, that taking the formal limit 
of the large-amplitude oscillations A;»Ao one obtains 
(rK)^r™"^(A)>r°. for the time-averaged Kondo tem- 
perature if r™°^<^U (see the insert in Fig. 2). In the 
opposite limit T Lfj{ti<t<t2)^U the system falls into a 
mixed valence regime where the Kondo temperature is 
the poorly-defined quantity. The result Q survives as 
long as the large amplitude limit holds provided that the 
condition (Tk) / Dq'^1 is still fulfilled. We conclude that 
the Kondo temperature considerably increases as com- 
pared to Tf^ when the shuttling particle approaches to 
one of the leads. The relative variation of the Kondo 
temperature oscillations at small shuttling amplitudes is 
given by 
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(5) 



where Xo^)^o/ \/h'i{DQ/T^)<$:_Xo is the effective tunnel- 
ing length which accounts for the Kondo rcnormaliza- 
tions. 
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FIG. 2: Differential conductance of a Kondo shuttle. Solid 
line denotes G for the shuttle, dashed line stands the for static 
island located in the symmetric (central) position, dot line 
stands for the static island located asymmetrically. Insert 
shows the time oscillations of Tk for small (dot line) and 
large (solid line) shuttling amplitudes. 

Let us discuss the temperature behavior of the differen- 
tial conductance G{T, Vdc -^ 0)=dI/dVdc- In the strong 
coupling Kondo limit T'^T^ 
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2e2 
h 



^5Tk\1 



T° 



X2 
^0 



.(6) 



The conductance does not reach the unitary limit Gjj = 
2e^//i due to the asymmetry in the respective couplings 
to the leads. If, however, the shuttling island was not 
centrally positioned when starting its motion, the effec- 
tive magnitude of the conductance can grow as compared 
to its value at the starting position (see Fig. 2). 



In the weak coupling regime T]J*"^^T<C£'o the zero 
bias anomaly (ZBA) in the tunneling conductance is 
given by 



G(r) = ^Ga( 



^TL{t)VR{t) 



I 



{Tat)+Tn{t)YMTlTK{t)W 



(7) 



Although the central position of the island is most favor- 
able for the Breit-Wigner (BW) resonance (Fi=rfl'), it 
corresponds to the minimal width of the Abrikosov-Suhl 
resonance. The turning points correspond to the maxi- 
mum of the Kondo temperature given by the equation (Q 
while the system is away from the BW resonance. These 
two competing effects lead to the effective enhancement 
of G at high temperatures: 



G{T) = Gl 



1 



I-2a2(T)sinh^(a;(i)/Ao) 



(8) 



where G%=Gu{?>T^'^)l{l&^ri^{T /T%)) is the conductance 
of the static island in the central position (Fig. 2), 
a^{T)=hi{Do/T)/hi{T/T^) and AT=Ao/a(r) is the tem- 
perature dependent tunneling length. Evaluating ^ for 
the small amplitude limit under the condition A^^CAq, 
we obtain 
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G{T) 






JT, 



K 



I 






Gl- 



TO ln(T/TO)- 



(9) 



Formally, the correction to the conductance 6Gk 
© must be compared with the regular term 
0(C/[ln^(r/T°)]). The latter, however has much 
smaller amplitude C-ln(ln(T/r]|.))<ln(Do/T]^). 
Thus, the Eq.@ describes the leading correction to 
conductance. 

In the limit T'^'^hfl<^T the differential conductance 
in the weak coupling regime is given by 



G 



■in' 



peak 



Gr 



I6"^[ln(ft/(rT0))P' 



(10) 



where h/r ^ Ml is determined by the dccoherence ef- 
fects associated with the non-adiabaticity of the mo- 
tion of the shuttling-particle and by the Q-factor of 
the NEM device. In general, the behavior of the dif- 
ferential conductance at low temperatures has a form 
Gj,,ak/Gu=F[[5TKlT%)-f{hn/T^j,)], where F{x) and 
f{y) are two universal functions, each of them depend- 
ing on one variable similar to [20j. In the large voltage 
limit eV^T^ the finite current transferred by the shut- 
tle and therefore the noise created by this current leads 
to h/r^eV destroying the Kondo effect. Wc will present 
more detailed discussion on the dccoherence effects in the 
" anti-adiabatic" Kondo regime elsewhere. 

Next we turn to the case of even A^ in the island (Fig. 
^b). In this case one may refer to the excited-state Kondo 



features [3, |2l| , where the KR tunnehng is possible only 
during the time intervals where 



Asrit) = S{t) ~ Je^t) < TK{t). 



(11) 



The level spacing 6{t)=e2{t)—ei{t) may reduce due to the 
tunneling-induced Friedel shift 



^^it) 



Y: \TiHt)\'Ref- 



pode 



(12) 



a=L,R 



.(2) 



.(1) 



provided Ta >Ta , which is usually the case |3|- This 
effect is maximal near the turning points of shuttle mo- 
tion. Similar 2-nd order tunnel processes results in the 
so called Haldane renormalization |22| of J ex- In close 
analogy with |2|, it is easy to see that the reduction of 
exchange gap in a dot with even occupation obeys the 
renormalization group flow equation 



dA5TM, = p„Vf|T(^)p-|TWp' 



(13) 



where rj=h\{DQ/D) is the scaling variable describing the 
reduction of the energy scale D of the band electrons 
in the leads. Additional contribution to this reduction 
originates from the mixture of the exited state with two 
electrons on the level £2 with the ground state singlet [3 ■ 
These effects are also maximal around the turning points 
of the shuttle trajectory. 

Thus, if the condition I|1HI is valid for the certain 
time intervals during the oscillation cycle (Fig.l.b), the 
Kondo tunneling is possible for a part of this cycle, where 
the shuttle is close to one of the leads. It should be 
emphasized that in this regime only the weak-coupling 
Kondo effect may be observed at T^T^, whereas at 
T^O the triplet state is quenched and the dot behaves 
as a zero spin nano-particle |2l|. The full scale Kondo 
effect may arise only if the variation of |Tq •^(i)p induces 
the crossover from a singlet to a triplet ground state of a 
shuttle. The singlet/triplet crossover induced by the vari- 
ation of gate voltages was observed on a static planar dot 
|23|. Since the adiabaticity condition /ifi<C T^ ~ ' is vio- 
lated close to the singlet /triplet transition, our approach 
does not apply to this regime. 

In excited-state Kondo regime Tk depends on 
the value of Ast, being scaled relative to the 



true S= 



-1 value 

,(S=1) 



(S=l) 



Tk 

(S=l) 



accordance 



with the law 



TK{t)lT)^-'> = [T)^-'>/AsT{t)\' with C<1 0. Thus, 
we conclude that the Kondo shuttling in case of even 
A'^ may be observed as a pulse ZBA in tunnel conduc- 
tion (Fig.l.c), which emerges in time intervals 5t, where 
the condition Hll|l is fulfilled. Assuming the linear de- 
pendence /S.st{.x) [23, we estimate these intervals as 
6t^6x/{ily/{x^)), where Sx is the distance from the turn- 
ing point at which the Kondo shuttling is possible |2J| . 

The Kondo shuttling differs from the Kondo effect in 
molecular conductor with the centrum of mass motion 



[l6j. In that case the shift of the centrum of mass 
is caused by single electron transport, and the non- 
adiahatic phonon-assisted processes interfere with the 
Kondo tunneling. 

In conclusion, wc have found that the Kondo shuttling 
in NEM-SET increases the Kondo temperature due to the 
asymmetry of coupling in the turning points compared to 
central position of the island. As a result, in the case of 
odd A^ the differential conductance is enhanced in the 
weak coupling regime and is suppressed in the strong 
coupling limit. In the case of even A^, Kondo tunneling 
exists only as a shuttling effect. 
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